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Abstract 

In order to adequately describe molecular rotation far from equilibrium, we have generalized the 
J-diffusion model by allowing the rotational relaxation rate to be angular momentum dependent. 
The calculated nonequilibrium rotational correlation functions (CFs) are shown to decay much 
slower than their equilibrium counterparts, and orientational CFs of hot molecules exhibit coher- 
ent behavior, which persists for several rotational periods. As distinct from the results of standard 
theories, rotational and orientational CFs are found to dependent strongly on the nonequilibrium 
preparation of the molecular ensemble. We predict the Arrhenius energy dependence of rota- 
tional relaxation times and violation of the Hubbard relations for orientational relaxation times. 
The standard and generalized J-diffusion models are shown to be almost indistinguishable under 
equilibrium conditions. Far from equilibrium, their predictions may differ dramatically. 
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I. INTRODUCTION 



There have been collected many evidences (due to polarization time-resolved 
experiments,-i^i2i>^i^'^'^i^'^'^»ii>i2i'i2i»i^'i^ computer simulations,— 'i^'i^'i^'^iSi and model 
studies^i^^iS^) which indicate that rotational and orientational relaxations of molecules in 
liquids under equilibrium and (highly) nonequilibrium conditions differ significantly. The 
studies of nonequilibrium rotational relaxation have recently culminated in papers.^"^ As 
has been demonstrated, relaxation of rotational energy of hot nonequilibrium photofrag- 
ments slows significantly with the increase of their initial rotational temperature and differs 
dramatically from relaxation of the equilibrium rotational energy correlation function, 
manifesting thereby breakdown of the linear response description. Furthermore, hot 
nonequilibrium molecules exhibit almost dissipation-free coherent rotation, which persists 
for several rotational periods. Molecular dynamic simulations have attributed this unusual 
behavior to a rapid rearrangement of the local liquid structure: a hot solute pushes away the 
neighboring solvent molecules and rotates almost freely until it looses enough energy.— 
Another way of thinking about that is in terms of the angular momentum dependent 
rotational friction: highly rotationally excited molecules experience lower friction than their 
thermally equilibrated counterparts.— 

These nonequilibrium relaxations are not reproduced by conventional theories of (equilib- 
rium) molecular rotation, which are based upon master equations with constant relaxation 
rates (the extended diffusion models,— '^'^'^ the Keilson-Storer model^S'^'^) or rotational 
Fokker-Planck equations with constant frictions.— i^i^i^ Physically, this is not surprising. If 
we consider relaxation under equilibrium conditions, then all the relevant energies are of the 
order of ksT {kB being the Boltzmann constant and T being the temperature of the bath). 
So, putting aside non-Markovian effects, we can always introduce a certain effective constant 
relaxation rate. If we consider relaxation under nonequilibrium conditions, when initial en- 
ergies are much higher than fc^T, such a description is no longer applicable, and the angular 
momentum dependence of relaxation rates must be explicitly taken into consideration (see 
Ref.— for more detailed argumentation and discussion). 

There exist two major and complimentary approaches to rotational relaxation, via the 
Fokker-Planck equations2^i^i^i2& and extended diffusion models.— i^*^"^ In our previous 
paper,^^ we have demonstrated how the angular-momentum dependent friction can be incor- 
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porated into the rotational Fokker-Planck equation. The present paper deals with a similar 
generalization of the J-diffusion model. First, we investigate how the angular momentum 
dependent rates manifest themselves in relaxations of initial nonequilibrium distributions 
and in behaviors of rotational and orientational correlation functions (CPs). Second, we de- 
velop a reliable tool for the description of nonequilibrium molecular rotation, which can be 
used i.e. for the interpretation of photofragment anisotropy decays.— ^'^'^'^^^'^^^^^'^^'^^^ 
Third, we demonstrate how the use of standard "equilibrium" approaches may lead to quan- 
titatively and even qualitatively wrong results, if they are apphed to describe nonequilibrium 
molecular rotation. 

Photodissociation ICN + hu CN + I (which was studied experimentally in the 
gas-phase^i^i^ and in the condensed phase^»i^, and also investigated via computer 
simulation3i^ii^'i^'2Si2ii2^) is chosen in the present article as a prototype process which pro- 
duces highly nonequilibrium hot photofragments. The results of the molecular dynamics 
simulations2ii2^ are used to obtain realistic values of the model parameters and to test our 
theoretical predictions. 

The structure of our paper is as follows. The generalized J-diffusion model, which ac- 
counts for the angular momentum dependence of the relaxation rate, is developed in Sec. 
2. A convenient explicit expression for the initial nonequilibrium angular momentum dis- 
tribution is derived in Section 3. Relaxation of this nonequilibrium distribution, as well as 
dynamics of the angular momentum and energy CPs, are investigated in Section 4. Orienta- 
tional relaxation is studied in Section 5. The Appendixes describe the methods of calculation 
of rotational CPs (A), their relaxation times (B), and orientational CPs (C). 

The reduced variables are used throughout the article: time, angular momentum and 
energy are measured in units of a/ // (kBT), y/IksT and /c_bT, respectively. Here fc^ is the 
Boltzmann constant, T is the temperature of the equilibrium bath, and / is the moment 
of inertia of the photofragment, so that = a/// {ksT) is the averaged period of its free 
rotation. Por CN at 300 K, Tr =0.2 ps. All the Laplace-transformed operators are denoted 
by tilde, viz. 




(1) 
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II. GENERALIZED J-DIFFUSION MODEL 



Restricting our consideration to linear molecules and neglecting the short-time non- 
Markovian effects, we start from the rotational master equation 

dtp{3, n, t) = -aLp(j, n,t)- J rfj'T(j'| j)p(J, ^,t) + j rfj'r(j| j')p(j', t). (2) 

Here p(J,r2,)f:) is the probability density function, J is the (two-dimensional) angular mo- 
mentum of the linear top in its molecular frame, Q, are the Euler angles which specify 
orientation of the molecular frame with respect to the laboratory one. The first term on the 
right-hand side of Eq. ([2]) is the free-rotor Liouville operator, which describes the angular 
momentum driven reorientation, L being the angular momentum operator in the molecular 
frame. The last two terms are responsible for the bath-induced relaxation. Since our molec- 
ular ensemble is isotropic, the kernel of the dissipation operator, T(J|J'), is il-independent. 
It must obey the detailed balance 

T(J|J')pb(J') =T(J'|J)pb(J), (3) 

which is responsible for bringing the system under study to the equilibrium rotational Boltz- 
mann distribution 

PB(J) = (2vr)-iexp{-jV2}. (4) 

Normalization of Eq. ([2]) , 

j d3dnp{3,n,t) = 1, (5) 

which insures conservation of probability, is accounted for automatically. 

By selecting a particular form of the relaxation kernel T(J'| J), one can recover many mod- 
els of rotational relaxation available in the literature: the extended diffusion models,— i^i^i^ 
the rotational Fokker-Planck equation,— i^*^*^ the Keilson-Storer model.— i^i^ All these 
models, however, yield j dJ'T{J'\3) = u = const, that is the collision rate is angular mo- 
mentum independent. For example, the most popular model of rotational relaxation, the 
J-diffusion model, corresponds to the choice 

T(J|J') = iypB{3). (6) 

We cannot straightforwardly generalize the J-diffusion model by replacing the rate u with 
its J-dependent counterpart i^{J), since this procedure would violate normalization (and 



therefore conservation) of the probabihty density (Eq. ([5])). To get a correct description, we 
assume that the relaxation kernel can be represented in the factorized Gaussian form 



T(J| J') = exp{-aJ72} exp{-6J'72}c, (7) 

a, b, and c being numerical constants to be specified. To satisfy the detailed balance ([3]), we 
must choose 

a = l + b. (8) 
Then the generalized J-diffusion relaxation kernel takes the form 

T(J| J') = u^-^ exp{-(l + 6) JV2} exp{-6J'V2}. (9) 
27r 

Here the rate u governs the dissipation strength and the constant 

(10) 

is chosen to insure that 

J dJ'T{3'\J) = vz{J), z{J) = exp{-6JV2}. (11) 

Thus, vzi^J) can be regarded as J-dependent collision frequency.— 

Eq. ([7]) looks as a natural generalization of the standard J-diffusion kernel ([6]), which still 
preserves its Gaussian and factorized form. If h is small (it is indeed small, see bellow), then 
the kernel ([7]) can be regarded as the first order correction to its J-diffusion counterpart ([6]). 
Furthermore, any relaxation kernel T(J'| J) can be expanded in series on Gaussians, 

T(J| J') = exp{-(l + 6) JV2} exp{-W'V2}, (12) 

b 

Uf, being certain expansion coefficients.^^ Thus, Eq. retains the first term in this expan- 
sion, and can further be refined by taking more terms, if necessary. The use of Eq. (JTj) can 
be justified a posteriori, since the predictions of the generalized J-diffusion model in Sees. 4 
and 5 sustain comparison with the results of molecular dynamics simulations performed in 

Plugging the kernel ([7]) into Eq. ([2]) we obtain our generalized J-diffusion master equation: 

dtp{J, n, t) = -aLp(j, n, t) - uz{j)p{3, n, t) 



+ z/(l + J)pb(J) J d3'z{J')p{J', n, t). (13) 

This is the equation which will be studied in the subsequent Sections. When b = {z = 1), 
it reduces to the standard J-diffusion model. 

Before embarking at particular calculations, it is useful to estimate plausible values of 
the parameters u and b. Since u describes the overall dissipation strength and is similar 
to its counterpart in the standard J-diffusion model, there are no intrinsic limitations on 
the value of this parameter. It is small (< 1) for rarefied gases and large (> 1) for liquids 
and solutions. As to the value of b, it is expected to be small (<^ 1), since the effects due 
to the J-dependence of relaxation rates do not manifest themselves at room temperatures 
under equilibrium conditions. On the other hand, if the initial nonequilibrium energy of the 
photofragment is much larger than fcsT, then the product b{E^/kBT) is not necessarily 
small. It is in this case the effects due to the J-dependence of relaxation rates become 
important. Note that the parameters u and b must be independent of initial conditions. 
However, as we shall see in Sees. 4 and 5, Eq. (fT3il predicts that nonequilibrium responses, 
i.e. the ensuing rotational and orientational CPs, depend strongly upon initial conditions. 



III. NONEQUILIBRIUM INITIAL DISTRIBUTION 

Normally, master equations like (JT3l) are used to calculate various rotational and orienta- 
tional CPs under equilibrium conditions, assuming that p(J,O,0) = {2h)^^pb{J)- We wish 
to study rotational relaxation under nonequilibrium initial conditions 

p(J,f2,0) = (27r)-Vne(J). (14) 

To get an explicit expression for pne{3), we adopt a model developed in Ref.~ Let us con- 
sider the photoreaction A + hu ^ B + products. We assume that the photofragmentation 
is instantaneous on the time scale of molecular rotation. This assumption holds true for 
most of sniall^'^'^'^'^'^'^ and polyatomic^ molecules, for which the photofragmentation 
time can be be estimated by few hundreds of femtoseconds.^ Having accepted the assump- 
tion about an instantaneous (impulsive) photodecomposition, we immediately conclude that 
dissociation produces a nonequilibrium distribution over the angular momenta of photoprod- 
ucts. Indeed, there exist two major sources of rotational excitation of fragments. These are 
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the parent molecule rotation and the applied torque. Therefore, the angular momenta of 
the parent (J^) and product (J^) molecules are connected through the formula^Ii^ 

JB,a= -^B,a-Ra/3(S)/^]3JA,/3 + ^a- (15) 

a=x,y,z 

The first term describes mapping of the parent molecule rotation into that of the product, 
and the angular momentum A originates from the impulsive torque arising due to the 
rupture of chemical bond(s) of the parent. The small Latin indexes label the Cartesian 
components of the corresponding vectors and tensors, lA,a and lB,a are the main moments 
of inertia of the species A and B, i?a/3(S) is the matrix of rotation from the frame of the 
main moments of inertia of the product to that of the parent, and S are the pertinent Euler 
angles. 

Let the parent molecule A be a planar asymmetric top (i.e., a triatomic molecule, Ia,x = 
lA,y + Ia,z) and the product B be a linear rotor {Ib,x = lB,y = Ib, Ib,z = 0). We then 
define the molecular frames of the parent (x^, za) and product (x^, ub, zb) molecules 
in such a way that the axes xaH^b are perpendicular to the plane of the parent molecule, 
and Zfi-axis coincides with that of the linear fragment (see Fig. 1). In this case 

f I \ 

cos(^ij) sin^ipB) 
y - sin(V^B) cos(V'b) J 

V'b = /-{yAllB) = ^{zaZb)- It is natural to surmise that the applied torque is perpendicular 
to the plane of the parent molecule, A||xb. Then, assuming that the parent molecules have 
a Boltzmann equilibrium distribution at the bath temperature T, we find that the linear 
photoproducts are distributed according tc^ 

Pne(J) = TT^ 

(hereafter, the subscript B, which denotes the angular momentum and moment of inertia 
of fragment B, is omitted and A is the magnitude of the impulsive angular momentum A). 
We have introduced the parameters^ 

I ' /A,x{cos2(V'B)//A,y + sin2(?/'B)//A,J 

Eq. (fTSi) is a desirable nonequilibrium distribution, which is employed in all the subse- 
quent calculations. The entire information about dissociation is contained in the parameters 



i?(H) 



(16) 



j t^^exp|-| ([Jcos^ + A]2 + eV2sin2^)} (17) 



A, ?7, and ^. The first of tliem is responsible for the apphed torque and can be considered as 
dynamical. The last two parameters control the parent-product rotational energy transfer. 
They are determined by the geometry of the dissociating molecule. Note that < ^ < 1. 
When ^ ^ 0, then the parent molecule is highly prolate {Ia,z ^ Ia,x, lA,y) and the linear 
fragment is "attached" perpendicularly to the parent molecule axis {jpB ^ /'^)- If ^ ~ 1, 
then the parent molecule is also considerably prolate, but the linear fragment is "attached" 
parallel to the parent molecule axis (^/'b ~ 0). If the parent molecule is a symmetric top 
{Ia,x = lA,y)) then ^ = 1. The quantities rj and ^ can be regarded as known for a particu- 
lar photoreaction, since the angle ipB and the main moments of inertia are fixed for specific 
molecules A and B. The value of A can be estimated through the consideration of the energy 
partitioning in the course of the photofragmentation.— 

Distribution f|T8|) is convenient for the further use, since it covers many different scenario 
of the photofragmentation. It reduces to the Boltzmann equilibrium distribution (jlj) if A = 
and 7] =^ = 1. On the other hand, we get the delta-distribution 

in the limit 1] oo. 

The quantum version of distribution (fTTI) is immediately written down after the "quan- 
tization" of the angular momentum:— 

oo 

plij) = ZQ'x{j), Zq = Y,x{j), (20) 

d9exp {-^ ([J cos^ + 5f + ef sin' o) } . 
Here j is the rotational quantum number, 

j = VjUTi), s = a/c, vq = C\ C^n/^ik^. (21) 

Eq. fl20l) allows us to simulate a rotational distribution of any width centered at any value 
of j. Several characteristic shapes of PneU) presented in Fig. 2. If ^ = 1, the distri- 
bution is symmetric. If ^ < 1 (for triatomics, this corresponds to dissociation from a bent 
configuration) the distribution becomes asymmetric, with a long-j bias. 

Let us now apply the above machinery to photodissociation ICN + /iz/ — > CN + I. This 
reaction proceeds via two channels, producing low energy and high energy photofragments. 

8 



The quantum distribution PneU) fo^^ nascent CN- fragments has been measured inM A com- 
parison of the "full line" distribution in Fig. 2 with that reported in™ reveals that the 
former reproduces qualitatively the experimental results for hot photofragments. For CN 
at 300 K, C = 0.135. Therefore, the "realistic values" of the dimensionless parameters 
V = Vq/C"^ = l-l) ^ = = 5.4 will be used in our illustrative calculations.— 



IV. ROTATIONAL RELAXATION 

In this section, we study evolution of the quantities which depend on the angular mo- 
mentum J but are independent of the Euler angles O. We thus can integrate Eq. ([T3|) over 
O and arrive at the reduced master equation 

dtp{J, t) = -uz{J)pi3, t) + u{l + b)z{J)pB{J) J dJ'ziJ')p{3', t). (22) 

This equation can be used, for example, to calculate the angular momentum CF 



the averaged rotational energy 



CAt) = (23) 

W /ne 



Csit) = <f'y-"<f^ (24) 



and the rotational energy CF 



^'^ ^ - (E''),^-(E)s(E)^ • * ' 

Here 

E = JV2 (26) 

and we use the notation = / dJpa{3)..., a = B, ne. The corresponding integral 

relaxation times are determined as 

POO 

Ta = / dtCait), a = J,E, S. (27) 



As is shown in Appendix A, the Laplace images of the probability density p(J, t), as well 
as of CFs (I23|) - (l25|) can be given in quadratures (Eqs. (lA2ip . (lASp . and (1A13|) ) for any 
initial nonequilibrium distribution p„e(J)- The inversion of the Laplace transforms into the 
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time domain can elementary be performed numerically. Furthermore, we derive elucidating 
analytical expressions in several important particular cases, which help us to grasp essential 
features of rotational relaxation under nonequilibrium conditions. 

Let us assume that the nonequilibrium distribution is given by Eq. flTIJl) . that is the 
photoproducts are produced with a fixed magnitude of the angular momentum A. Such a 
situation corresponds to the procedure of the (microcanonical) preparation of hot molecules 
in simulations reported We thus define the nonequilibrium photofragment energy 

= AV2. (28) 

We assume, in addition, that b is small {b <^ rj, t]^'^) but the product bE/^ can take any 
value. Then Eq. (]A20p predicts that the CFs (123|) - (!25|) are all the same and exponential, 

Cait) = exp{-uz{A)t}, a = J, E, S. (29) 

The J-dependent rates (fTTj) are seen to slow rotational relaxation. The corresponding re- 
laxation times (|27|) read 

^« = = - exp < > . (30 

The last term in this expression is written in dimensional units, to emphasize that the 
generalized J-diffusion model predicts the Arrhenius energy dependence of the rotational 
relaxation times. Eq. fl5Ul) corroborates the finding^i2^'2^ that there exist characteristic 
"nonequilibrium energy" E/^ and "nonequilibrium temperature" Ta = E/^/ks, 

E^r^kBT/b, Ta~T/6 (31) 

starting from which the J-dependence of rotational relaxation becomes important. If the 
characteristic nonequilibrium energy E/^ is much larger than ksT, then the relaxation times 
for hot photofragment s can substantially be longer than their equilibrium counterparts = 
1/z/, despite b is small. The standard J-diffusion model, on the contrary, predicts that any 
rotational CP flAll) decays exponentially with the rate z/, irrespective of particular forms of 
A{J), -B(J), and Pne{J) (see Eq. (IA19P ). It is thus inadequate under highly nonequilibrium 
conditions, when 6i?A > 1- 

The Arrhenius energy dependence of the rotational relaxation times correlates with the 
results of Tao and Stratt,— who simulated CFs Cs{t) (12^ for different initial nonequilib- 
rium temperatures Ta. They obtained the following values: (a) ts = 5.29 for Ta = 1917 
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K, (b) Ts = 7.73 for Ta = 2395 K, (c) ts = 11.25 for Ta = 2875 K (see also Ref.-^). We 
can take any two sets of ts and Ta to calculate u and b by using Eq. ( l30i) . Doing so for 
three different pairings of the sets, we have got remarkably consistent results: u = 0.87, 
b = 0.095 for sets a&b, u = 0.84, b = 0.094 for sets b&c, u = 0.86, b = 0.095 for sets c&a. 
If we recall that Eqs. fl29|) and fl30l) are approximate and valid in the limit of b <^ 1, rj ^ 1, 
the agreement becomes even more encouraging. Thus the quantity z{A) is close to 1 (0.91) 
at equilibrium (Ta = 120 K), 0.22 at Ta = 1917 K , 0.15 at Ta = 2395 K, and 0.11 at 
Ta = 2875. So, rotational relaxation at Ta = 2875 K is nine times slower than that at the 
equilibrium temperature. Summarizing, the value of 6 = 0.1 is a "realistic value" for the 
parameter which is responsible for the J-dependent relaxation. 

It is timely to present an exact expression for the angular momentum integral relaxation 
time, which is valid for any values of the parameters of the generalized J-diffusion model. 
According to the formulas derived in Appendix B, 



{v{b) is defined via Eq. (1B8|) ). If we assume that b r], rjC,"^, then Eq. (l32i) reduces to (130|) 
for any r], ^, and A. This lends an additional support to the validity of Eqs. f l30l) . 

On the other hand, Eq. (l32l) demonstrates that the J-dependence of rotational relaxation 
rates works both ways: if A is small and r] < 1, then nonequilibrium integral relaxation 
times can be smaller then their equilibrium counterparts. This is illustrated by Fig. 3, 
which shows the reduced angular momentum relaxation time tj/tj'^ (full line) and rotational 
energy relaxation times ts/t^'^ (dashed line) vs. parameter rj. The nonequilibrium relaxation 
times are calculated via Eqs. (]B4|) - (1B9I) . and their equilibrium counterparts are given by 
Eqs. (IB 121) . The observed decrease of the relaxation times at A ^ 1 and r] < 1 is not so 
pronounced as their exponential increase at A ^ 1. 

Fig. 4 compares rotational CFs ( I23l) -( l25l) calculated within the standard (6 = 0) and 
generalized (6 = 0.1) J-diffusion models. The CFs have been computed as explained in 
Appendix A. The realistic values of the model parameters (r7 = l.l,^ = l,A = 5.4, and 
z/ = 1) have been used. The standard J-diffusion predicts all the CFs to coincide with 
exp{— z/t}. According to the generalized J-diffusion model, the CFs relax much slowly, as 
expected. For the model parameters chosen, Cs(t) and Cj(t) are almost indistinguishable,— 




(32) 
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while CE{t) decays slower than the former two CFs. 

If we assume that the initial distribution is given by delta-function f|T9|) and b <^ r], t]^"^, 
then the general Eqs. (]A2ip and (lA22p simplify to yield the following evolution of the 
probability density function: 

p(J, t) = ^ ^^^~ exp{-uz{A)t} + pb{3){1 - exp{-uz{A)t}). (33) 

In an accord with what has been observed in,— the rotational distribution fl55]) is clearly 
bimodal. At every time moment, p(J, t) is a mixture of the initial nonequilibrium distribution 
(fT9|) and the equilibrium distribution (j4]). The standard J-diffusion model predicts the same 
formula but with -2(A) = 1. Thus, the J-dependent rate effectively slows the rotational 
relaxation and increases the lifetime of the initial nonequilibrium contribution. The higher 
is the value of bE/^, the stronger is the effect. For example, if we take the dimensionless values 
of the parameters u = 0.84, b = 0.094 obtained above, we can calculate the decay rates at 
Ta = 2875 K. Namely, the generalized J-diffusion model predicts z/z(A) = 0.3 ps^^, while the 
J-diffusion model yields u = 2.6 ps^^. The generalized J-diffusion rate explains why p{J,t) 
simulated for CN fragments at Ta = 2875 K shows both the nonequilibrium and equilibrium 
contributions for more than 4 ps.^^ The standard J-diffusion model, which predicts the decay 
rate to be 9 times higher, is thus absolutely inadequate far from equilibrium. 



V. ORIENTATIONAL RELAXATION 

Orientational CP of the rank j is defined through the Wigner D-functions^ as follows: 

G'it) =< D^{n{t))D^{n{0))-^ > . (34) 

Its Laplace image, G^{s), can be calculated as explained in Appendix C, Eq. (lC8p . If we take 
the initial delta- function distribution f|T9|) . then Eq. flC8P can considerably be simplified. At 
short times (t < (z/z(A))^^), it can be inverted into the time domain to yield 

G^{t) ^ exp{~uz{A)t} ^Ojo + ^ 5Z "j'^ cos(A;At) j , (35) 

the numerical coefficients ajk are given by Eq. flC5p . The standard J-diffusion predicts 
the same formula but with 2; (A) = 1. Thus, the oscillations in orientational CFs of hot 
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photofragments, which have been measured in the gas phase^i^^ and simulated in the con- 
densed phase,— are entirely determined by the initial nonequilibrium distribution ([T9l) . 
The period of these oscillations is uniquely determined by the value of A. 

Fig. 5 compares orientational CFs simulated for hot CN fragments injected at the tem- 
perature Ta = 2875 K into the heat bath of argon atoms at T = 120 K (ReisM'^) and those 
calculated within the generalized J-diffusion model. The delta-function initial distribution 
f|T9|) has been used in our calculations, since it corresponds to the procedure of the prepa- 
ration of the ensemble of photofragments in the simulations. The simulated orientational 
CFs look qualitatively very similar to those described by Eq. ( l35i) . Quantitatively, they 
can be fitted by Eq. ([35D relatively well, but the so-obtained decay rate of the second rank 
orientational CF turns out to be 1.4 higher than that of the first rank CF. The orientational 
CFs which are exactly calculated within the generalized J-diffusion model are seen to re- 
produce the simulated CFs very well. This is remarkable, since we did not fit the simulated 
CFs: we used the parameters of the generalized J-diffusion model obtained in the previous 
Section through the comparison of the simulated and theoretical energy relaxation times. 
Namely, we took b = 0.094, u = 0.84, and A = 6.75.— This lends an additional support to 
self-consistency and predictive strength of the generalized J-diffusion model. 

Eq. (135|) makes it clear that the persistence of the oscillatory behavior for orientational 
CFs is much higher in the generalized J-diffusion model, due to the A-induced decrease 
of the relaxation rate, z/z(A). This is vividly illustrated by Fig. 6, in which presented 
are the first and second rank orientational CFs calculated within the generalized (b = 0.1) 
and standard (6 = 0) J-diffusion models for the "realistic values" of the model parameters 
(?7 = 1.1, ^ = 1, A = 5.4, and) and u = 7 . The generalized J-diffusion model predicts the 
highly oscillatory orientational CFs. In exactly the same situation, the standard J-diffusion 
model predicts overdamped and slowly decaying CFs, which have nothing in common with 
the "true" CFs. Thus, the use of the standard (equilibrium) models of rotational relaxation 
beyond their domain of validity may lead to completely wrong predictions. 

We can use Eq. ( ]C8[) to derive a simple expression for orientational CFs in the case of 
strong dissipation, z/ ^ 1. As expected, they are described by the diffusion formula 

G^{t)=exp{-Dj{j + l)t} (36) 

with the diffusion coefficient D equaled to the equilibrium angular momentum integral re- 
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laxation time, rj^. The latter is given by Eq. (lB12p . so that 



The orientational relaxation times are calculated as 



oo 



= / dtG^{t) = G^{0). (38) 
Jo 

It is popular to plot these quantities vs. the angular momentum relaxation times or rates.— 
That is why Figs. 7 display the graphs r^^ (a) and Tq (b) vs. u calculated within the standard 
and generalized J-diffusion models. 

To get a better idea about the behavior of these curves, it is insightful to obtain explicit 
expressions for in case of weak and strong dissipation. If <^ 1, then Eq. (ICSP predicts 
that Tj^ is inversely proportional to the rate u: 

= ^ f ^^^_exp(;^j + ii±^^ . (39) 

Here Cj = ajo/(2j + 1) and ajo are explicitly given by Eq. (IC5]).-^ The standard J- 
diffusion formula^^i^ is recovered at 6 = 0. Again, if the product bEj\ is ~ 1 or larger, 
then nonequilibrium initial conditions manifest themselves in a considerable increase of 
in comparison with the standard J-diffusion predictions ai u < 1 (Figs. 7b). 

When u > 1, the situation changes and the standard J-diffusion model overestimates the 
actual values of (Figs. 7). If ^ 1 then the diffusion formula holds and, therefore, 

Arf = — ^— . (40) 

nj ^^^^^^ V ) 

This expression can be termed as the generalized Hubbard relation. The standard Hubbard 
relation, t^tj = [j{j + l)]~^,i^i^i^ is seen to be significantly off. This discrepancy is re- 
markable, since it embodies the breakdown of the linear response theory. Indeed, the onset 
of rotational diffusion occurs after the molecular angular momenta have been thermolized 
according to the bath- induced Boltzmann distribution (jlj). Thus the rotational diffusion 
equation (15^ is a legitimate description a.t t > tj and the rotational diffusion coefficient 
is determined by the equilibrium rotational fluctuations and therefore by the equilibrium 
value of the angular momentum relaxation time, rj^. As has been detailed in Section 4, 
if the system starts far from equilibrium, then its angular momentum relaxation time tj 
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can differ dramatically from its equilibrium counterpart, Tj*^. It is in this case we expect 
the breakdown of the linear response Green-Kubo-type formulas for transport coefficients, 
which identify the rotational diffusion coefficient D with the angular momentum relaxation 
time Tj. 

VI. CONCLUSION 

An adequate description of rotational and translational relaxation in liquids under 
nonequilibrium conditions (or within an interval of characteristic energies which highly 
exceeds that of the bath thermal energies) cannot be carried out in terms of conven- 
tional Langevin and Fokker-Planck equations with constant frictions or by master equa- 
tions with constant relaxation rates: one has to take into account (linear and/or angu- 
lar) velocity dependence of friction.— i^iiS^iSSii^i^i^i^M The present paper deals with study- 
ing nonequilibrium rotational and orientational relaxation. We have generalized the stan- 
dard J-diffusion model by allowing its dissipation rate to be angular momentum depen- 
dent and calculated various nonequilibrium rotational and orientational CFs. The reaction 
\CN + hv —>■ CN + I has been selected as a prototype process which produces highly nonequi- 
librium hot photofragments. We have used the results of computer simulations^^i^^ to obtain 
realistic values of the model parameters and to test our theoretical predictions. 

We have found that nonequilibrium rotational relaxation rates assume the form 
u exp {—bE^/{kBT)}, where u is the equilibrium rate, is the initial "nonequilibrium 
energy" , T is the equilibrium temperature of the heat bath, and b is the dimensionless small 
parameter fixed for a system under study. Accordingly, the relaxation times have the Ar- 
rhenius energy dependence, exp {bE^/ (ksT)}. So, there exist characteristic "nonequi- 
librium energy" E^ ~ ksT/b and "nonequilibrium temperature" Ta ^ T/b, starting from 
which the J-dependence of relaxation rate becomes important and induces a significant slow- 
ing of rotational relaxation. In agreement with simulations,—*^ relaxation of the rotational 
probability density is shown to be a bimodal time-dependent mixture of the initial nonequi- 
librium and final equilibrium distributions. The slowing of rotational relaxation induces a 
long lifetime of the initial nonequilibrium contribution, which manifests itself in pronounced 
coherent effects. 

The slowing of rotational relaxation causes qualitative changes in molecular reorientation. 



15 



In agreement with simulations,— hot nonequihbrium molecules exhibit slightly perturbed 
coherent rotation, which persists for several rotational periods. A similar observation has 
recently been made in the context of molecular excitation by strong femtosecond pulses: 
nonequihbrium molecular wave packets can be formed in such a way as to slow their subse- 
quent rotational relaxation.'^ 

We have demonstrated that orientational CFs in the overdamped limit are described 
by the diffusion equation with the diffusion coefficient equaled to the equilibrium angular 
momentum integral relaxation time, ry'. Thus the orientational relaxation times obey 
the generalized Hubbard relation, r^rj*^ = + 1)]^^. Since the nonequihbrium angular 
momentum relaxation time tj can differ substantially from its equilibrium counterpart Tj'^, 
the standard Hubbard relation,— i^*^ in which rj*^ is replaced by tj, can be significantly off. 

It should be emphasized that all rotational and orientational CFs and their relaxation 
times dependent explicitly on the nonequihbrium preparation of the molecular ensemble. 
This effect is unreproducible within the standard (equilibrium) rotational models, which 
predict that relaxation rates are independent of initial conditions and are given by the linear 
response theory. Thus, the Arrhenius forms of the rotational relaxation times, slowing down 
of rotational relaxation, and violation of the Hubbard relations are all manifestations of the 
breakdown of the linear response theory far from equilibrium. 

In practical terms, the standard and generalized J-diffusion models are almost indistin- 
guishable under equilibrium conditions. Under nonequihbrium conditions, their predictions 
differ dramatically. The differences are not only quantitative but, not infrequently, qualita- 
tive. The message is thus as follows: the friction and/or the relaxation rate must be taken 
J-dependent in order to adequately describe rotational relaxation far from equilibrium. 

Note, finally, that our theory can straightforwardly be generalized to symmetric and 
asymmetric top molecules. In this latter case, the relaxation rate can be taken rotational 
energy dependent, rather than angular momentum dependent. The method of solution of 
the corresponding kinetic equations remains absolutely the same, and the explicit formulas 
for various CFs can be written down after a straightforward generalization of the results 
presented in the Appendixes. Our theory is also readily extendable to quantum case, by 
replacing the integrations over the angular momentum J with summations over the rota- 
tional quantum number j, and switching from the classical orientational CFs fICSP to their 
quantum mechanical counterparts (see, e.g., Ref.— ). Such a generalization might be use- 
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ful for describing rotational coherences and time-dependent alignments in nonequilibrium 
dissipative systems, like those studied in Ref.— 



APPENDIX A: CALCULATION OF ROTATIONAL CFS 

Let A{J) and B{J) be arbitrary functions of the angular momentum. Then the CF 

CABit) = {A{3{t))B{3))^^ (Al) 
can be evaluated through the rotational kinetic equation fl22p as follows: 

CABit) = j dJA{3)p{3,t), p(J,0)=p„e(J)5(J). (A2) 
After applying the Laplace transform ([1]) to Eq. (!22l) with initial condition (lA2p . we obtain 

Here 

a{s) = j d3z{J)p{J,s). (A4) 

Multiplying Eq. (1A2|) by z{J) and integrating over J, we obtain an algebraic equation for 
Inserting its solution into Eq. flA3l) . multiplying the so obtained expression by A{J) 
and integrating over J, we get 

Cab{s) = <!>ne{AB, s) + u{l + h) ' > ^^\ ' (A5) 

Here we have introduced the functionals 

$„(y,t) = y"rfJp„(J)exp{-z/^(J)t}F(J), a = B,ne (A6) 
and their Laplace transforms 

$,(y,g)= / (iJp,(J) a = B,ne (A7) 

J s -\- VZyJ ) 

which are defined for any function V(J). 

In principle, Eq. f lASp delivers the desirable expression for CF Cab{s) which, after the 
numerical inversion of the Laplace transform, yields the CF in the time domain, CAB{t)- 
The long-time limit of this CF is determined as follows: 

CAB{t^c^)^{A)j,{B)^^. (A8) 
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Here the averages are defined as 



= v^a{Yz, 0) = J d3p,{3)Y{3), a = B, ne. (A9) 

If both {A)^ and (-B)„g are nonzero, then CF (lAip possesses a stationary long-time asymp- 
tote flASp . This is so for the rotational energy CF, for example. Therefore, Cab{s) — »• 
(A)^ (-B)„e when s — > 0. This is undesirable for doing numerics. It is more convenient to 
subtract this constant contribution and redefine the CF ( lAlj) as follows: 

= {^B)„-(A>b(B>„.- 
Evidently, this CF is normalized to unity (Xa_b(0) = 1) and does not have any stationary 
asymptote {XAB{t) when t — oo). 

Taking the Laplace transform of Eq. (lAlOl) , we get 



Eq. (lAlip possesses yet an undesirable property: its numerator is a difference of two terms 



~ Of course, these two terms cancel each other when s ^ but it is convenient for 
numerical implementations to explicitly extract this singular contribution out of Cab{s). To 
this end, we make use of the identity 

^a{Yz,s) = ^[{Y)^-s^aiY,s)], a = B,ne, (A12) 

which is immediately derived from Eq. flA7l) . Applying this formula several times, we can 
rewrite Eq. (lAllI) in the following equivalent form 



Here 

5 (s) = {A)^ {B)^^ s) + s^b{A, s)^UB. s) - {A) J, s) - $^(^4, s). (AM) 



Eqs. (IA13P and flA14p are our final singularity-free formulas, which are used for the 



numerical inversion of the Laplace transforms. The explicit formulas for CFs ( l23l) - (!25l) can 
be obtained by inserting the corresponding functions A{3) and -B(J) into Eqs. (]A13I) and 
(IA14p . For example, the angular momentum CF (1231) corresponds to A[3) = B{J) = J. 
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Due to the isotropy of rotational space, ^a{.Jz"'{J), s) = 0, so that the general expression 
( 1A13P simplifies considerably: 

CAs) = (A15) 

There is no such a simplification for CFs ([23]) and (1251) . and all the terms (lA14p contribute 
into Eq. flAlSp . We do not give the corresponding expressions explicitly, since they add 
nothing profound. 

CF (lAlOp can be used to define an important quantity, the integral relaxation time: 

POD ^ 

TAB= dtXABit) = XAsiO). (A16) 

Jo 

Since u^b{z,0) = 1, Eqs. flAlSp and flA14D predict that 



Here 



vl/(0) = {A)^ (B)^^ $^(1, 0) - {A)^ 0) - 0). (A18) 

If we assume that the relaxation rate is constant {b = 0, z{J) = 1) then the J-diffusion 

model is recovered. In this case, evidently, 

(Y) ~ (Y) 
$,(r,s) = ^^, $,(F,0) = ^, a = B,ne. 
s + u u 

Then, Eqs. flM3|) and flXTTD simplify to 

Xab{s) = {v + s)~\ XabW = exp{-z/t}, tab = {^^T^ (A19) 

irrespective of particular forms of A(J), -B(J), and p„e(J)- Thus all rotational CFs in the 
J-diffusion model decay exponentially with the rate z/, and their integral relaxation times 
are all equaled to the inverse value of this rate. 

If the relaxation rate is J-dependent (6 7^ 0), the situation is much more complicated. 
However, the formulas simplify dramatically if the nonequilibrium distribution is given by 
the delta-function flTI?]) . If we assume, additionally, that 6^1 but the product hE^^ can 
take any value, we can write then 



s + ' s + uz{/\) 
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After the insertion of these expressions into Eqs. (1A13P and (lA14p . we obtain 



Xab{s)= , XAB(t) = exp{-z/z(A)t}, tab = (A20) 

for any A(J) and -B(J). 

If we wish to follow how the initial nonequilibrium distribution p„e(J) transforms in time 
into the equilibrium one, pb(J), we must solve kinetic equation (!22l) with the initial condition 
p(J, 0) = Pne(J)- The solution is given by a slightly modified version of Eqs. ( ]A3[) and ( lASj) : 

P^^^ = 7T^) ' ^^^^^ 

a,{s) = ^id^ . (A22) 

APPENDIX B: CALCULATION OF ROTATIONAL INTEGRAL RELAXATION 
TIMES 

If we employ the nonequilibrium initial distribution fll7p . then the quantities 0) 
(Eq. (]A7|) ) can be evaluated analytically for any CP of interest in the present paper 
{Y = 1, E, E"^). This means that the integral relaxation times tab can also be calculated 
analytically. To this end, convenient is to introduce the generating function 

T{b)= j JdJpne{J)eMKJl + Jy)l'^} = 

^JdJ.J dJ^exp{-7]{J,,-Af/2}exp{-7]{^Jyf/2}exp{b{4 + J^)/2}. (Bl) 
Evidently, 

~ 1 d^T(h) ~ ~ 

$.e(i^;^0) = -^^, <l>B(i?^0) = $„e(i^^,0)|5=,=l,A=0; (B2) 

(E^)^ = z/$,(E^,0)|fe=o, a = B,ne. (B3) 
Eq. (IBip is elementary evaluated to yield 

_ i ybA-" \ 

Differentiating this expression with respect to b, we obtain: 

$„e(l,0) = ^T(6), (B5) 
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<!>ne{E,0) = ^T{b)v{b), (B6) 



<l>ne{E\0) = ^T{b) {v\b)+w{b)} . (B7) 



Here 



These expressions can be substituted into Eq. (1A17I) to calculate various integral relaxation 
times.— In general, the so-obtained expressions are quite cumbersome and are not presented 
here. However, they are useful for obtaining simple elucidating formulas in several particular 
cases. If 6 -C r], r]^'^ and A 1, then 

T(6) ^ exp 1^ j , v{b) ^ — , w{b) ^ — . (BIO) 

If the initial distribution (IT7I) reduces to the equilibrium one (?7 = ^ = 1, A = 0), we get 

m = v{b) = w{b) = (Bii) 
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and^- 

z/(l-6)2' z/(l-6)= 
It has not escaped our notice that the integral relaxation times tab ( 1A17I) diverge for 
b > T], as predicted by Eqs. ( ]B4l) -( !B9l) . Although such situation is far beyond the expected 
domain of validity of our model (6 ^ r^), there is nothing pathological in the occurrence of 
the "phase transition" at 6 = r/ and all CPs are well behaved for b > rj. This means simply 
that the J-dependence of relaxation rates slows rotational CPs so dramatically that their 
integral relaxation times do not exist. 

APPENDIX C: CALCULATION OF ORIENTATIONAL CFS 



Having inserted the definition of orientational CP ( IMI) into Eq. ( 1131) . we obtain the 
equation 

dtG^{J,t) = -t3VG'{3,t) - iyz{J)G^{J,t) 
21 



+ u{l + b)z{J)pB{3) j d3'z{J')G^{J',t). (CI) 

Here are the matrix elements of the angular momentum operators over the D- 
functions:— 

{Li)ki ± liL^ki = 4,iTi{(j ±1)UtI + 1)V^\ {Li)ki = I5ku -3 < kj < j. (C2) 
Eq. (ICip must be solved with the initial condition 

G^(J,t = 0) = p„e(J). (C3) 

Evidently, 



G^{t) = J dJG^{J,t). (C4) 

Orientational CPs are calculated very similarly to rotational CPs (see Appendix A). Pirst 
we introduce the free linear rotor orientational CP— 

j 2 

F^(J,t) = + 2^ajfcCos{A;Jt}, ajk = (dl^{^)^ , (C5) 

k=l 

^kmiP) being the reduced Wigner function.— We further define the functional 

Q{(Y,t) = j d3pai3)exp{-uz{J)t}Y{3)F^{J,t), a = B,ne (C6) 
and its Laplace transform 

QiiY, s) = j dJp,(J)r(J)F^( J, s + uz{J)) = (C7) 

k=-j 

Then, closely following the derivation of Eq. flASp . we obtain the following expression for 
the Laplace transform of the orientational CP: 

&is) = QUI, s) + Ki + &)7^^7r^|M7^- ^^8) 

1 - iy{l + bjQ'siz^ s) 

Since orientational CPs of dissipative molecules, as distinct from their bath-free counterparts 
(1C5I) and rotational CPs (Appendix A), do not possess stationary asymptotes, this formula 
is well behaved in the limit s ^ and can be used for the numerical inversion of the Laplace 
images into the time domain. 
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FIG. 1: Sketch of the dissociating triatomic molecule. 
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FIG. 2: Rotational distributions for hot photofragments calculated for rjq = 0.02, ^ = 1, (5 = 40 
(full line), r]Q = 0.02, ^ = 0.2, S = 40 (dotted line). The equilibrium Boltzmann distribution for 
CN at 300 K is depicted by dashed line (r?Q = 0.018, ^ = 1,6 = 0). 
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FIG. 3: Reduced angular momentum relaxation time tj/tj'^ (full line) and rotational energy re- 
laxation time Ts/Tg^ (dashed line) vs. parameter r] calculated in the generalized J-diffusion model 
(5 = 0.1) for = 1 and A = 0. The dotted line corresponds to the situation when the reduced 
value of the relaxation times equals 1. 
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FIG. 4: Angular momentum CFs Cj{t) (full lines) and rotational energy CF CE{t) (dashed line) 
calculated for the "realistic values" of the model parameters {rj = 1.1, ^ = 1, A = 5.4) and u = 1. 

The lower line corresponds to the standard J-diffusion model (6 = 0) and the upper lines correspond 
to the generalized J-diffusion model (6 = 0.1). 
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FIG. 5: Orientational CFs of the first (a) and second (b) rank. Full lines depict the results of 
molecular dynamics simulations of CN fragments injected at the temperature Ta = 2875 K into 
the heat bath of argon atoms at T = 120 K (Refs.— and dotted lines show the results of 
the calculations within the generalized J-diffusion model (6 = 0.094, v = 0.84) with the initial 
^-distribution ([ED centered at A = 6.75. For CN at 120 K, =0.32 ps. 
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FIG. 6: Orientational CFs of the first rank (full lines) and second rank (dashed lines) calculated 
for the "realistic values" of the model parameters (77 = 1.1, ^ = 1, A = 5.4) and u = 7. The upper 

lines correspond to the standard J-diffusion model (b = 0) and the lower lines correspond to the 
generalized J-diff'usion model (6 = 0.1). 
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FIG. 7: Orientational relaxation times Tq of the first (a) and second (b) rank vs. the collision 
rate u calculated for the "realistic values" of the model parameters (r/ = 1.1, ^ = 1, A = 5.4). The 
dashed lines correspond to the Tq and modernized Hubbard asymptote l/ljij + 1)tj^] calculated 
within the standard J-difFusion model (6 = 0). The full lines correspond to their generalized J- 
difFusion model (6 = 0.1) counterparts. The standard Hubbard asymptotes, + 1)tj], are 

depicted by dotted lines. 
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